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Evaluation of the importance of spin-orbit couplings in the nonadiabatic quantum 
dynamics with quantum fidelity and with its efficient "on-the-fly" ab initio 

semiclassical approximation 

Tomas Zimmermann and Jin VanicekQ 

Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingenierie Chimiques, 
Ecole Polytechnique Federale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland 

We propose to measure the importance of spin-orbit couplings (SOCs) in the nonadiabatic molec- 
ular quantum dynamics rigorously with quantum fidelity. To make the criterion practical, quan- 
tum fidelity is estimated efficiently with the multiple-surface dephasing representation (MSDR). 
The MSDR is a semiclassical method that includes nuclear quantum effects through interference 
of mixed quantum-classical trajectories without the need for the Hessian of potential energy sur- 
faces. Two variants of the MSDR are studied, in which the nuclei are propagated either with the 
fewest-switches surface hopping or with the locally mean field dynamics. The fidelity criterion and 
MSDR are first tested on one-dimensional model systems amenable to numerically exact quantum 
dynamics. Then, the MSDR is combined with "on-the-fly" computed electronic structure to measure 
the importance of SOCs and nonadiabatic couplings (NACs) in the photoisomerization dynamics of 
CH2NH2 + considering 20 electronic states and in the collision of F + H2 considering six electronic 
states. 



I. INTRODUCTION 

Nonadiabatic couplings (NACs) originating in the 
Born-Oppenheimer separation of motion of electrons and 
nuclei often play an important role in the molecular 
dynamics^ In photochemistry, NACs are responsible for 
radiationless decay of the electronically excited states by 
a process called internal conversion^ In chemical reac- 
tion dynamics, they may affect branching ratios of prod- 
uct channels or allow for reactivity of otherwise non- 
reactive states, when the transition from a nonreactive 
potential energy surface (PES) to a reactive PES is in- 
duced by the motion of nuclei 4^ In some cases, transi- 
tions between PESs may be induced by another type of 
couplings called spin-orbit couplings (SOCs), which cou- 
ple the electronic spin with angular momentum. 6 Rela- 
tivistic in origin, SOCs become increasingly important 
in molecules containing heavy elements. Nevertheless, 
SOCs often play a role also in molecules composed of light 
elements, e.g., various chromophores or DNA bases. 7 
In photochemistry, SOCs cause radiationless transitions 
between PESs of different multiplicity called intersys- 
tem crossings, which are responsible for the phospho- 
rescence or emergence of triplet intermediates in pho- 
tochemical reactions^— Similarly to NACs, SOCs may 
affect chemical reactions even when no photoexcitation 
is involved. 11 ' 12 

Even though intersystem crossing is typically slower 
than internal conversion in molecules composed of light 
elements, in some cases both effects may occur on a com- 
parable time scaled - — Therefore, there is a need for 
a rigorous criterion of the importance of SOCs in ac- 
curate quantum nonadiabatic simulations. Knowing in 
advance which PESs and which couplings are impor- 
tant for given initial conditions may speed up a simu- 
lation and avoid costly calculations of additional PESs. 
Below, we propose and test a criterion of importance 



of SOCs which is based on the quantum fidelity^ 7 - In 
a similar manner, fidelity was already used to measure 
nondiabaticityi 8 - or nonadiabaticit y 19 ' 20 of quantum dy- 
namics. However, direct evaluation of fidelity requires 
knowledge of the quantum state evolved with the com- 
plete Hamiltonian including SOCs and consequently is 
impractical as a guide for quantum simulations. Fortu- 
nately, there exists a fast semiclassical method to com- 
pute quantum fidelity, called the multiple- surface dephas- 
ing representation (MSDR) , 18 ' 20 Evaluation of the im- 
portance of SOCs with the MSDR is the main focus of 
this work. 

Originally, the MSDR was introduced to measure ei- 
ther the nonadiabaticity or nondiabaticity of the quan- 
tum molecular dynamics. 1 ' Roughly speaking, "nondi- 
abaticity" is the difference between the diabatic quantum 
dynamics (i.e., quantum dynamics in which the diabatic 
coordinate couplings between diabatic electronic surfaces 
are neglected) and the fully coupled quantum dynamics. 
Similarly, "nonadiabaticity" is the difference between the 
adiabatic quantum dynamics (i.e., quantum dynamics in 
which the nonadiabatic momentum couplings between 
adiabatic electronic surfaces are neglected) and the fully 
coupled quantum dynamics. In both cases the "differ- 
ence" between the two types of dynamics is measured by 
the decay of the overlap of the time-dependent molecular 
wave functions evolved using the two types of dynamics. 
In the same manner, the importance of SOCs may be 
evaluated either in: 1) the basis which is diabatic with 
respect to SOCs, meaning that the Hamiltonian matrix 
expressed in this basis does not contain any spin-orbit- 
related couplings that depend on the nuclear momen- 
tum, or 2) in the basis which is adiabatic with respect 
to SOCs, meaning that the only spin-orbit-related cou- 
plings present in the Hamiltonian depend on the nuclear 
momentum. The basis which is diabatic with respect to 
SOCs is usually more suitable when SOCs represent only 
a small perturbation, whereas the basis which is adia- 
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batic with respect to SOCs is more appropriate in case 
of strong coupling. Independently of the adiabaticity or 
diabaticity with respect to SOCs, the basis may be adia- 
batic or diabatic with respect to couplings originating in 
the Born-Oppenheimer separation of motion of electrons 
and nuclei. Below, we use exclusively the basis which is 
diabatic with respect to SOCs. When analyzing the im- 
portance of SOCs, we always evaluate the importance of 
the full spin-orbit Hamiltonian which may include diago- 
nal terms as well. Nevertheless, the fidelity criterion and 
MSDR may be used also in the fully adiabatic basis of 
spin-orbit states, which is employed, for example, in the 
mixed quantum-classical code SHARCj^i It is important 
to note that spin-orbit-related couplings have different 
meaning in the two basis sets. In the diabatic basis with 
respect to SOCs, when the full spin-orbit Hamiltonian 
is considered as a perturbation, SOCs reflect faithfully 
the importance of the spin-orbit interaction in the sys- 
tem. In the adiabatic basis with respect to SOCs, on 
the other hand, the spin-orbit interaction modifies the 
PESs themselves and the lack of importance of the offdi- 
agonal spin-orbit-related couplings does not necessarily 
mean that the spin-orbit interaction is unimportant. 

The MSDR is a generalization to nonadiabatic dynam- 
ics of the dephasing representation (DR))22r— derived 
for dynamics on a single surface using the Van Vleck 
propagator or the linearization of the path integral^ In 
the single-surface setting, the DR is closely related to 
the semiclassical perturbation approximation of Smith, 
Hubbard, and Mille r 26 ' 27 and to the phase averaging of 
Mukamelj2£ Its applications include evaluations of the 
stability of quantum dynamic a 23 ' 24 ' 29- — and calculations 
of electronic spectra ^ 8 ' 32 " — 

While the MSDR is not a method for general dynamics, 
it can be used to compute all quantities which can be ex- 
pressed in terms of quantum fidelity amplitude. Weaker 
generality and semiclassical nature give the MSDR some 
advantages in comparison to general methods for nona- 
diabatic dynamics. 

The main advantage of the MSDR compared to 
wave packet methods is that the number of trajecto- 
ries needed for convergence of MSDR does not scale ex- 
ponentially with the number of degrees of freedom. 38 
Therefore, the MSDR may be applied to problems 
with dimensionality beyond the scope of even the 
most advanced methods for nonadiabatic quantum dy- 
namics such as the multi-configuration time-dependent 
Hartree methods^ - — Similarly to other semiclassical 
methods J 3 -^ the MSDR includes some quantum effects 
on the nuclear motion. The advantage of MSDR in com- 
parison to most other semiclassical approaches is that 
the MSDR does not require the Hessian of the potential 
energy, which is often the most expensive part of semi- 
classical calculations (see, e.g., Ref. 51). In terms of 
the computational cost per trajectory, the MSDR thus 
falls into the category of methods in which nuclei are 
treated classically such as surface hopping methods^ - — 
and other methods based on the mixed quantum-classical 



Liouville equation^ - — or methods obtained by lineariza- 
tion of the path integral expression for the quantum 
propagator. 64 ' 65 (Costs of current implementations of the 
MSDR are essentially the costs of a mean field or surface 
hopping dynamics. Nevertheless, in principle, the MSDR 
may be combined with a propagation scheme based di- 
rectly on the mixed quantum-classical Liouville equation, 
which is typically more expensive but also more accu- 
rate.) 

The outline of the paper is as follows: SectionlHl starts 
by introducing the fidelity criterion of importance of 
SOCs and by defining the spin-orbit Hamiltonian. This 
is followed by a brief derivation of the MSDR, and the 
description of the propagation scheme and the algorithm. 
The section ends with the computational details. In Sec- 
tion lllll the fidelity criterion and MSDR are tested us- 
ing one-dimensional model systems which allow for nu- 
merically exact quantum solution. Subsequently, the 
MSDR is combined with "on-the-fly" computed ab ini- 
tio electronic structure and applied to evaluate the im- 
portance of SOCs and NACs in the photoisomerization 
of CH 2 NH 2 + and in the collision of F + H 2 . SectionEl 
concludes the paper. 

II. THEORY 

A. Fidelity as a measure of importance of 
spin-orbit coupling terms 

Following our work on the nondiabaticity 18 and 
nonadiabaticity 20 of the molecular dynamics we base the 
quantitative criterion of importance of SOCs on quantum 
fidelity Fqm between molecular quantum states propa- 
gated with and without SOCs. More precisely, 

F Q M(t) = \f Q M(t)\ 2 = \(r(tW(t))\ 2 , (1) 

where /qm is the quantum fidelity amplitude, |^°(i)) = 
e -iH°t/fi |^,(o)) is the quantum state of the molecule 
evolved using the nonadiabatic Hamiltonian H° which 
does not include SOCs of interest and \ip € (t)) — 
e -iWt/h |^(o)) is the quantum state evolved using the 
fully coupled nonadiabatic Hamiltonian H £ = H° + eV. 
In the last expression, V = H so contains SOCs of in- 
terest and e controls the extent of perturbation. (Bold 
face denotes n x n matrices acting on the Hilbert space 
spanned by n electronic states, hat * denotes nuclear op- 
erators. In general, the superscript or e designates 
the Hamiltonian with which the object was propagated.) 
When F QM « 1, \ip°(t)) is close to \^{t)) and SOCs do 
not influence the dynamics significantly. On the other 
hand, when Fqm •C 1, SOCs are important and should 
be taken into account in an accurate quantum calcula- 
tion. The advantage of the fidelity criterion is that, in 
addition to population transfer between PESs (which is 
a standard dynamical measure of the importance of cou- 
plings) , fidelity can detect subtle effects on the dynamics 
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caused by the displacement and interference on a single 
PES (see Ref.H3 for details). 



B. Spin-orbit Hamiltonian 

In principle, any spin-orbit coupling Hamiltonian H so 
may be used. In this work, model potentials are used in 
one-dimensional systems. In case of the photoisomeriza- 
tion of CH2NH2 + and collision of F + H2, elements of 
H so are computed with the Breit-Pauli Hamiltonian 66 
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where qjj — (qj — Qj) is the difference of position vec- 
tors of electron j and nucleus J, qjk = (qj — qk) is the 
difference of position vectors of of electron j and k, Pj is 
the momentum of electron j, and Sj is the spin of elec- 
tron j. Exceptionally, since H BP is not expressed in any 
basis yet, in Eq. © (and only there) hats denote both 
electronic and nuclear operators. 



C. MSDR 

Here we will only summarize the theory of MSDR; the 
full derivation may be found in Ref.HTj. The derivation 
starts with a quantum fidelity amplitude expression gen- 
eralized to the Hilbert space given by the tensor product 



/ QM (*) = Tr ( 
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where p mlt is the density operator of the initial state. 
Expressing /qm in the interaction picture and partially 
Wigner transforming 6 ^ the resulting equation over nu- 
clear degrees of freedom yields an alternative exact ex- 
pression 



/qm (f) 



h- D Tr e 



j2D „init 
d Ap w 



(X) ■ {re' i ^o^ l ( t ') dt ' /h 

(4) 

where X denotes the point (Q, P) in the 2D-dimensional 
nuclear phase space, Tr e is the trace over electronic de- 
grees of freedom, and V 1 (t) is the perturbation in the 
interaction picture given by H°. The first approximation 
consists in replacing the Wigner transform of a product 
of operators in the Taylor expansion of the time-ordered 
exponential in Eq. (U) by the product of Wigner trans- 
forms of these operators. Recognizing the resulting sum 
as an Taylor expansion of another exponential, this ap- 
proximation can be expressed succinctly as 



To evaluate this expression, the time evolution of 
V^y (X, t) has to be known. The second approxima- 
tion involves replacing the exact evolution by a mixed 
quantum-classical (MQC) propagation scheme described 
below, leading to the final expression for the MSDR of 
fidelity amplitude, 

/msdr(<) = (re""'" v^ >MQC (x, t >'M 
= h- D Tr e J d 2D X P ^ t (X)-Te' l,! ^ v ^,MQc(x,t')dt'/h^ 

(6) 

D. Propagation scheme 

The MQC equation 2 ^^! 2 . for the evolution of the 
density matrix is given by 

<9Pw,MQC « r„ 1 
qI - ~^ L^W, Pw.MQCj 

+ o ({ H W,Pw,MQc} - {Pw,MQC> H w}) , (7) 

where the explicit dependence of p w mqc on time and 
on the nuclear phase-space coordinate X was omitted 
for clarity. The propagation equation for M qc (X, t) , 
which is the last ingredient needed in our method, dif- 
fers from Eq. ([7J only by the sign of the time deriva- 
tive. Therefore Eq. © together with Eq. ([7]) define the 
MSDR. Several numerical approaches exist that solve 
Eq. in terms of "classical" trajectories X(t). How- 
ever, since trajectory-based methods for solving Eq. ([7]) 
are still relatively complicated, the MSDR is in practice 
implemented using one of two schemes which further ap- 
proximate Eq. ((7j). The common feature of the two ap- 
proximations is that all elements of p w (X, t) are prop- 
agated using the same PES (which may, nevertheless, 
differ for different trajectories). The first approach 2 ^ ap- 
proximates Eq. ([Jj as 
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where (A) e = Tr e (p e ■ A) is a partial average of A over 
the electronic subspace. This approach was called the 
locally mean field dynamics (LMFD) in Ref.l20. where it 
was derived simply by invoking the locally mean field ap- 
proximation. As pointed out in Ref.HTJ, the LMFD turns 
out to be nothing else than the dynamics of a swarm of 
trajectories, each of which is propagated with the Ehren- 
fest dynamics. The second approach employs the physi- 
cally motivated fewest-switches surface hopping (FSSH) 
scheme, 52 where the matrix elements of the density oper- 
ator p w fssh are_[in the adiabatic basis) computed using 
Eq. (II) of Ref.il. 
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E. Algorithm 

The details of the general implementation of the 
MSDR algorithm are given in Ref.|2y. In this work, we 
consider only initial states for which the "conditional" 
electronic density matrix is pure for all X and hence the 
initial density matrix can be written as the tensor prod- 
uct 

p™* (X) = p init (X) c init (X) ® c init (A) f (9) 

where p init {X) := Tr e fffi (A) and c init {X) is the ini- 
tial electronic wave function for nuclei located at X. In 
practice, JVtraj phase space points are sampled from a 
scalar nuclear density p lnlt (X) and a vector c lmt (X) is 
attributed to each of the generated phase space points. 
In the case of LMFD, this determines the initial condi- 
tion completely. In the case of FSSH, one also needs 
to select the initial surface randomly and separately for 
each trajectory according to the following prescription: 
For a trajectory starting at X, the probability for its ini- 
tial surface to be surface j is given by |c™ lt; (X)| . The 
trajectories are then propagated with [because, ul- 
timately, we propagate M q C (X, t) in the interaction 
picture given by H^] by using either the LMFD or FSSH 
dynamics. It is advantageous that for states ©, when 
the LMFD or FSSH dynamics is used, Eq. ^ simplifies 
to the weighted phase space average 

/msdrW = (c°(X,rf -c e (A,i)) . (10) 

In Eq. (|10p , the wave function c° (X, t) is obtained au- 
tomatically from the LMFD or FSSH dynamics. Analo- 
gously, c £ (X, t) solves the Schrodinger equation 

^^ = -1hw(x°W)-cW) (ii) 

for a single discrete electronic degree of freedom with a 
time-dependent Hamiltonian H^-(t) := H^ v (X°(t)) (i.e., 
in the Lagrangian reference frame given by H^) and with 
the initial condition c £ (A,0) = c init (X). In Eq. (dTJ), 
X°(t) is the phase space point resulting from the evolu- 
tion of the initial phase space point X for time t with 
H^-, using either the LMFD or FSSH dynamics. 



F. Computational details 

All quantum calculations were performed using the 
second-order split-operator algorithm. 73 The LMFD or 
FSSH dynamics were done using the second-order sym- 
plectic Verlet integrator! 7 ^ The Schrodinger equation 
for the discrete "electronic" degree of freedom was 
solved using the unitary propagator U(A, t,t + At) = 
e -«H w (x)At/fi._ js^ initio calculations were performed us- 
ing Molpro 2010.1Z& and Columbus 5.9.2^ 



III. RESULTS 

A. Comparison of the MSDR with the numerically 
exact quantum dynamics 

In order to test the utility of quantum fidelity as a mea- 
sure of the importance of SOCs and the performance of 
the MSDR in comparison to the numerically exact quan- 
tum dynamics we use three one-dimensional model sys- 
tems. Results are shown in Fig.ffl 

The PESs and couplings of model A are shown in the 
right part of Fig.[T](a). Both Hamiltonians H° and H c are 
expressed in the diabatic basis with respect to both SOCs 
and couplings originating in the Born-Oppenheimer sepa- 
ration of motion of electrons and nuclei. Both Hamiltoni- 
ans contain three PESs, of which the lower two are identi- 
cal to the PESs of Tully's single avoided crossing model. 52 
The third PES is flat with constant energy E — 0.15 a. u. 
In both Hamiltonians, the lower two PESs are coupled 
by the same coupling term V\% as in the original single 
avoided crossing model. Additionally, in H e (but not in 
H°), the highest PES is coupled to the lower two PESs 
with Via = ^3*1 = ^23 = ^3*2 = (1 + *) [Cexp (-DQ 2 )] , 
where C — 0.005 and D = 1.0. A more general, complex 
form was chosen to emulate spin-orbit coupling terms, 
which may also be complex-valued. The initial state is 
a Gaussian wave packet (GWP) located on the lowest 
energy PES with the mean kinetic energy To = 0.025 
a.u. As can be seen in Fig.[T](a), quantum fidelity -Fqm 
decays substantially (by more than 50%), signifying the 
importance of V13 and V23 in the dynamics. The decay of 
Fqm is accurately reproduced by the MSDR using both 
the LMFD and FSSH dynamics. On the other hand, the 
decay of the survival probability P1+2.QM = 1 — is 
very small and the final probability J^qm of finding the 
system on the third PES is less than 1%. This demon- 
strates that the survival probability -Pi+2,qm may be a 
poor measure of the importance of couplings between 
PESs. More refined criteria, such as the fidelity criterion, 
are needed to measure the influence of SOCs. (Note that 
in this specific case, strong effects of couplings of the first 
two PESs to the third PES may be inferred from consid- 
ering the probabilities Pi and P2 separately. However, 
for that it is necessary to run the dynamics twice - once 
with couplings to the third surface and once without the 
couplings. To use the fidelity criterion approximated by 
the MSDR only one simulation is sufficient.) 

Model B shown in Fig.Q](b) is based on the two-surface 
Tully's double avoided crossing models It is also ex- 
pressed in the diabatic basis with respect to both SOCs 
and couplings originating in the Born-Oppenheimer sep- 
aration of motion of electrons and nuclei. The third 
additional surface is described by the equation V33 = 
FQ 2 + G, where F = -2 • 10~ 5 and G = 0.06. The cou- 
pling term V23 = V£ 2 = ii exp (- JQ 2 ) (with I = 0.003 

and J = 0.001) is present only in H e . As can be seen in 
the right part of Fig.[T](b), the coupling term V23 is rela- 
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tively weak but widely spread, resembling closely SOCs 
typically seen in molecular dynamics. The initial state 
is a GWP located on the second lowest-energy PES with 
mean kinetic energy To = 0.368 a.u. Again, the decay 
of Fqm is closely followed by both implementations of 
the MSDR. This time, even P1+2.QM follows Fqm rela- 
tively well. Still, the approximate MSDR method gives 
a slightly more accurate picture of the influence of SOCs 
on the dynamics. 

Model C [Fig.Q](c)], inspired by Tully's extended cou- 
pling model^ demonstrates that the fidelity criterion 
detects the importance of SOCs even in cases where no 
significant transition of the probability density between 
surfaces occurs. The PESs and couplings are given by 

V u (Q) = -A 
V 22 (Q>0) = A + B cxp (-CQ) 
V 22 (Q<0) = A + B[2- cxp (CQ)} 
V 12 [Q > 0) = B [2 - exp {-CQ)] 
V12 {Q <Q) = B exp (CQ) , 

where A = 10, B = 0.1, and C = 0.9. The real-valued 
coupling term V12 is present in H e but not in H°. (Note 
that when the coupling term is purely imaginary, the re- 
sults are exactly identical.) As can be seen in Fig.Q](c), 
for a low energy wavepacket the wavefunctions propa- 
gated with H° and H e are very different. Indeed, fidelity 
decays towards zero. On the other hand, the survival 
probability Pi never decays by more than 6 • 10" °. Even 
though the model is slightly artificial (in that the PESs 
are considerably coupled, despite a large energy differ- 
ence between them), it clearly demonstrates the plausi- 
bility of situations in which the dynamical importance of 
SOCs would be undetectable if measured with the sur- 
vival probability as a criterion. Finally, Fig.Q](c) shows 
that the MSDR still closely follows the quantum result. 
[Since H° is uncoupled, both the LMFD and FSSH dy- 
namics are equivalent. Therefore, only the FSSH result 
is shown in Fig.Q](c).] 

B. On-the-fly ab initio photoisomerization 
dynamics of CH 2 NH 2 + 

As the first on-the-fly ab initio application of the 
MSDR, we examined the importance of SOCs and NACs 
in the quantum dynamics of the second excited sin- 
glet state of CH 2 NH 2 + using the SA5,5-CASSCF(6,4)/6- 
31G* method. The time step of the FSSH dynamics 
was equal to ~ 0.2 fs. In total, five lowest-energy sin- 
glet states and 15 lowest-energy triplet states (five states 
for each Ms = 1,0,-1) were considered in the simula- 
tion. To evaluate the nonadiabaticity, the unperturbed 
Hamiltonian H° did not contain any couplings between 
surfaces and the perturbed Hamitonian H e contained 
NACs between singlet states. To evaluate the impor- 
tance of SOCs, H° contained only NACs between singlet 
states while H c contained NACs between singlet states 



and SOCs between all 20 considered surfaces. The ini- 
tial state was the vibrational ground state of the ground 
electronic PES computed in the harmonic approxima- 
tion. (The harmonic approximation was used only to 
compute the initial state, and not for the propagation 
itself.) This wave packet was placed on the second ex- 
cited PES (in the basis diabatic with respect to SOCs 
and adiabatic with respect to couplings originating in the 
Born-Oppenheimer separation of motion of electrons and 
nuclei). Immediately after excitation, the system quickly 
decayed to the first excited state and subsequently to 
the ground state. Three main pathways were observed 
during the first 100 fs, in agreement with experimental 
observations and previous calculations £Zr2& 1) photoiso- 
merization leading to the hot ground state CH2NH2 + 
2) bi-pyramidalisation leading to dissociation into CH2 + 
and NH2 and 3) release of H2 . 

Fast decay of fidelity -Fnac, which can be seen in Fig. [51 
demonstrates the vast importance of NACs in the quan- 
tum dynamics of photoexcited CH2NH2 + . On the other 
hand, SOCs are not important in this dynamics as signi- 
fied by the very slow decay of fidelity -Fsoc ■ 



C. On-the-fly ab initio collision of F + H2 

The collision of F + H2 represents a system where 
SOCs are known to be important at least at low 
energieSi 12 ' 80 In order to check the ability of the MSDR 
to detect the importance of SOCs, we have applied it 
together with on-the-fly computed electronic structure 
using the SA3-CASSCF(9,10)/6-31+G* method. The 
(9,10) active space was chosen because it is known to 
produce qualitatively correct smooth PESs for any orien- 
tation of H2 axis4i The time step of the FSSH dynamics 
was equal to ~ 0.2 fs. Six lowest-energy doublet states 
were considered in the simulation (three states with 
M s = 1/2 and three states with Ms = -1/2). The ini- 
tial state was placed either on the electronic ground state 
or on the first excited state with Ms = 1/2 in the basis 
diabatic with respect to SOCs and adiabatic with respect 
to couplings originating in the Born-Oppenheimer sepa- 
ration of motion of electrons and nuclei. (The degeneracy 
of p states of F was removed by H2 at a distance of ~ 6 
A.) The molecule of H2 was in the vibrational and rota- 
tional ground state with uniformly random orientation of 
the molecular axis. Translational degrees of freedom of 
F and H2 were sampled from GWPs with the mean total 
kinetic energy of the collision Ek = 0.7 kcal/mol. At such 
a low energy, both considered collisions were nonreactive. 

As can be seen in Fig.|3l fidelity due to SOCs decays 
significantly in both cases. Relatively fast decay of fi- 
delity in the case of electronically excited initial state 
signifies importance of SOCs on the collision dynamics. 
In the case of the ground state collision, fidelity decays 
more slowly and SOCs may probably be neglected in less 
accurate simulations. Nevertheless, the decay is sufficient 
to justify inclusion of SOCs in accurate quantitative sim- 
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Figure 1. Importance of SOCs in the nonadiabatic dynamics: comparison with the numerically exact quantum dynamics on 
model systems. Left panels compare the quantum fidelity (Pqm) with the MSDR (Pmsdr) and with the quantum survival 
probability (Pqm)- Right panels show the corresponding diabatic PESs and couplings Vjj(Q) as well as initial [V>(0)] and final 
[i/>°(i/) and ip e (tf)] wave functions evolved with H° and H e , respectively. In all three models, the initial state was the GWP 
Tp(q) = -^j^ exp [— (Q — Qo) 2 /2a 2 + Po(Q — Qo)/ti] with the mass equal to 2000 a.u. (a) Coupling terms V13 = V23, containing 
both real and imaginary parts, result in a very low decay of Pqm which does not correspond to the substantial decay of 
Pqm- Only the lowest-energy PES was occupied initially by a GWP with Qo = — 15 a.u., Pn = 10 a.u., and a = 2.83 a.u. 

(b) Extensive but relatively weak purely imaginary spin-orbit coupling term causes considerable decay of both Pqm and Pqm- 
Only the second lowest-energy PES was occupied initially by a GWP with Qo — — 15 a.u., Pn = 38.36 a.u., and a — 0.74 a.u. 

(c) The coupling term V12 causes decay of Pqm towards zero without affecting Pqm- Only the lowest-energy PES was occupied 
initially by a GWP with Qo = —20 a.u., Pn = 0.2 a.u., and a = 3.41 a.u. 
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SOCs and NACs in the collision of F + H 2 [SA3-CASSCF(9,1 OJ/6-31 +G 
1 



0.8 



0.6 



0.4 



0.2 




/ 'grnd_SOC 
f exc_SOC 

F grnd_NAC 
F exc_NAC 



20 



40 



60 



80 



100 



f(fs) 




Figure 2. Importance of the SOCs and NACs in the dynamics 
starting on the second excited electronic state of Cfi2NH2 + 
evaluated with the MSDR combined with the FSSH dynamics. 
Top: Fast decay of FirAC signifies the importance of NACs, 
whereas a very slow decay of Jsoc shows that SOCs may 
be safely neglected. Bottom: Swarm of 72 trajectories rep- 
resenting the wave packet during the photoisomerization dy- 
namics of Cfi2NH2 + immediately after excitation and ~ 16 fs 
after excitation. Configurations on the first, second, and third 
lowest-energy excited states are depicted with yellow, red, 
and pink hydrogen atoms, respectively. Configurations on 
the ground state have grey hydrogens. 



Figure 3. Importance of SOCs and NACs in the collision of 
F + H 2 evaluated with the MSDR combined with the FSSH 
dynamics. Top: The decays of i ? grn d_soc and F exc _soc sig- 
nify the importance of SOCs, especially in the collision ini- 
tiated on the first excited electronic state. NACs for both 
excited and ground state dynamics are as important as SOCs 
for the ground state dynamics. Bottom: Swarm of 64 trajec- 
tories representing the wave packet (initially located on the 
first excited PES) during the collision of F + H2 at t = fs 
and t ~ 48 fs. Configurations on the first and second lowest- 
energy excited doublet states are depicted with yellow and red 
hydrogen atoms, respectively. Configurations on the ground 
state have grey hydrogens. 



illations. On the other hand, the decay of fidelity due 
to NACs is similar for both states. For the ground-state 
dynamics, NACs are comparable in importance to SOCs, 
whereas for the excited state dynamics, SOCs are clearly 
much more important. 



IV. CONCLUSIONS 

We have demonstrated that the fidelity criterion may 
detect disturbances of the quantum dynamics due to 
SOCs, which may not be found if cruder criteria, such as 
the extent of the population transfer between PESs, are 
used. This observation is in accordance with our previ- 
ous findings on the nondiabaticity and nonadiabaticity of 
quantum dynamics ; 18,20 It should be noted that fidelity is 
a nonspecific criterion and more specific measures may be 
constructed when needed in order to attribute the effect 
of couplings to either the electronic or nuclear dynamics. 
This can be done, for example, by combination of the 
fidelity criterion and surface population criterion. How- 
ever, the separation of the nuclear and electronic effects is 
not always possible. Fidelity is a rigorous single measure 



that can take into account both effects simultaneously. 

We found that the MSDR approximation of quantum 
fidelity remains accurate when the perturbation is caused 
by SOCs. Based on our previous experience, the MSDR 
may fail, especially in cases where the underlying mixed 
quantum-classical dynamics does not approximate the 
quantum dynamics reasonably well. This can happen, 
e.g., when tunneling is important in the dynamics of the 
unperturbed Hamiltonian H°. In that case, any imple- 
mentation of the MSDR, which is based on the mixed 
quantum-classical dynamics, may not work. Other fail- 
ures are specific to current implementations of MSDR, 
where both the LMFD and FSSH dynamics may be in- 
accurate due to the fact that all matrix elements of the 
density matrix are attached to the same trajectory and 
evolve on the same PES. A typical situation where ap- 
proximations of this kind fail occurs when regions of 
strong coupling are encountered repeatedly and dynam- 
ics on PESs between these encounters differ substantially. 
Surprisingly, in many cases the MSDR stays accurate de- 
spite the failure of the underlying dynamics. This is due 
to the cancellation of errors between the dynamics of H° 



s 



and the (not explicitly performed) dynamics of H e . 

Applications of the MSDR to the photoisomerization 
dynamics of CH2NH2 + and to the collision of F + H2 
demonstrated that our method may be used as a tool to 
decide which PESs and couplings are important in the 
quantum dynamics, without the need to run quantum 
dynamics itself. In contrast to the quantum dynamics, 
the MSDR does not scale exponentially with the number 
of degrees of freedom, it may be used on the fly and 
it does not require costly scans of PESs or any other 



substantial prior knowledge of a system. (Note that there 
exist several methods of quantum dynamics which may 
also be used on the fly.— ~— Nevertheless, they are not 
currently widely used.) Finally, the MSDR is simple and 
its current formulations may be very easily implemented 
into any code for FSSH or Ehrenfest dynamics. 
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